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Abstract A study of directional solidification of a 
weak binary alloy (specifically, Bi - 1 at% Sn) based 
on the fixed grid single domain approach is being 
undertaken. The enthalpy method is used to solve for 
the temperature field over the computational domain 
including both the solid and liquid phases; latent heat 
evolution is treated with the aid of an effective 
specific heat coefficient. A source term accounting 
for the release of solute into the liquid during 
solidification has been incorporated into the solute 
transport equation. The vorticity-stream function 
formulation is used to describe thermo-solutal 
convection in the liquid region. 

In this paper* we present a numerical simulation of 
g-jitter. A background gravity of 1 jig has been 
assumed, and new results for the effects of periodic 
disturbances over a range of amplitudes and 
frequencies on solute field and segregation have been 
presented. 

INTRODUCTION 

The investigation of solidification and melting 
processes in low gravity conditions is of great 
practical importance for crystal growth techniques. 
The quality of single crystals grown from the melt 
strongly depends on growth morphology and macro- 
segregation in solidified ingots caused by convective 
heat and mass transfer effects. A low gravity 
environment produces conditions in which convection 
is decreased to a level at which crystal growth is 
largely diffusion controlled. Residual accelerations in 
orbiting space vehicles are of the order of one to 
several hundred jig (where 1 pg =9.81xl0' 6 ms’ 2 ). For 
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this reason, much effort has been expended in recent 
years in performing crystal growth experiments in the 
microgravity environment of a spacecraft in earth 
orbit. Such effects as compositional and kinetic 
supercooling, and the influence of convection on 
compositional distribution in the melt, have been 
investigated under microgravity. 

However, gravity in an orbiting space vehicle may not 
be steady in either magnitude or direction. 
Perturbations to the anticipated steady microgravity 
environment may arise from, for example, crew 
actions, the operation of machinery and thruster rocket 
firings. Such perturbations are known as g-jitter. 

The effects of gravity perturbations on composition 
distribution in Bridgman crystal growth configuration 
have been investigated numerically in a number of 
papers. Alexander, Ouazzani and Rosenberger 1 
investigated the effects of steady and impulse residual 
acceleration on dopant distribution in Bridgman- 
Stockbarger crystal growth with different gravity 
vector orientations. It was found that lateral non- 
uniformity in composition is very sensitive to the 
orientation of the steady component of the residual 
gravity vector. 

Transient and periodic accelerations have been 
considered in Alexander, Amiroudine, Ouazzani, and 
Rosenberger 2 . It was found that the largest 
compositional nonuniformities occur for disturbances 
with amplitudes above 10' 6 g and frequencies below 
10' 2 Hz. At higher frequencies, larger acceleration 
amplitudes are required to obtain significant 
nonuniformities. Numerical results for the effect of 
g-jitter on the average interface concentration during 
Bridgman crystal growth in space are presented in 
Garandet, Corre, Favier and Alexander 3 . In all of 
these works a pseudo-steady-state model was adopted 
with the constraint that the interface is planar. This is 
a simplification of the true unsteady solidification 
process. 



In actual growth situations the solid-liquid interface 
can be non-planar due to thermal and mass transfer 
conditions and also due to morphological instability 
factors. Pseudo-steady-state model neglects transient 
effects such as changes in the velocity, temperature 
and concentration with time due to changes in the 
length of melt caused by translating the ampoule. 
Investigation of solute redistribution during the initial 
transients becomes crucial for an alloy with a low 
partition coefficient solidifying at low rates because 
the steady state is difficult to reach in laboratory 
experiments. 

In this work we investigate effects of periodic gravity 
perturbations on segregation and solute distribution 
during transient directional solidification of Bi-l%Sn 
alloy in a Bridgman furnace. A background gravity of 
I jig, which corresponds to a typical spacecraft 
environment, is considered. 


^ + V.(vc)=OV ! C ( 4 ) 

where t, p, p, c p , X and D are respectively the time, 
density, viscosity, specific heat and thermal 
conductivity of the alloy and the diffusivity of the 

solute; £ y/ t T t V and C are respectively the 
vorticity, stream function, temperature, velocity vector 
and solute concentration; g is the magnitude of the 
gravitational acceleration, and g is the unit vector in 
the direction of gravity. The density in the buoyancy 
term of equation (1) is assumed to be a linear function 
of temperature and solute concentration: 

P = P*[1-Pt(T-T m )+P c (C-C m )] (5) 


MATHEMATICAL FORMULATION 

We consider a Bridgman furnace in which a moving 
temperature profile consisting of a cold zone (T c ), a 
nominally adiabatic zone and a hot zone (T h ) is 
imposed on the boundary of the ampoule. This 
boundary temperature profile is translated with a 
constant pulling velocity as a result of the furnace 
movement causing the solid/liquid interface to move 
along the ampoule. The material in the ampoule is 
thus divided into two sub-regions: solid and liquid. 

Although the ampoule is three-dimensional, a two- 
dimensional model is used. This simplification is 
valid because, under the microgravity conditions being 
considered, convection is very weak and the 
solidification process remains largely diffusion- 
controlled. Newtonian and laminar flow is assumed in 
the liquid phase, and the Boussinesq approximation 
has been used, in which the liquid density is assumed 
to be constant except in the buoyancy term of the 
equation of motion. 


where j} T and fi c are the thermal and solutal 
expansion coefficients, 


and 


Pr=~ 


1 dp 
Pr 


p c =±*e. 

p,3C 


( 6 ) 

(7) 


p R is the reference density (the density of pure Bi at its 
melting temperature T R ) and C R is the Sn 
concentration in pure Bi, namely zero. 

The gravitational acceleration is taken as: 

g( O = go + A sin( 2na)t ) (8) 

where A is the amplitude of the acceleration, cuis the 
frequency and g 0 is the steady component of the 
acceleration. 


The governing time dependant equations describing 
mass, momentum, heat and solute transport in the 
vorticity-stream function formulation are: 

p ^r + v P x £|s| + ^v J f (i) 

VV = “C (2) 


Enthalpy method 

To model the process of directional solidification we 
have chosen the enthalpy method 4 which avoids 
explicit tracking of the solid/liquid interface. 

Latent heat evolution during phase change is 
incorporated in the energy equation using the 
following definition of enthalpy. For each phase 0, 
enthalpy is defined as 


pc p 
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where L is latent heat and /, is the local liquid volume 
fraction. 

For isothermal phase change the liquid fraction is 
determined by the melting temperature T m : 


for T >T m /,= 1 
for T<T m // = 0 


( 10 ) 


With the assumption that specific heat c p<> is constant 
in each phase, (9) can be written as 


where 2AT is a temperature interval chosen to 
represent the range over which phase change occurs in 
the (i, j) control volume. 

Based on the calculated values of liquid fraction at 
each mesh point the computational domain was 
subdivided into sub-regions of solid and liquid phases. 
In the solid, the vorticity, stream function and 
velocities were set to zero. In the liquid, they are 
calculated from the stream function defined as: 


V = Vxv? 


(16) 


h = c p ^T + f t L = h sens + f t L , 


(ID 


Solute transport with phase change 


Here h sens is the sensible heat, and the subscripts l and 
s refer to the liquid and solid phases. 

Using the apparent heat capacity method 5 , an effective 
specific heat can be defined by 


c-m=— =c, +& 
dT p * oT 


02 ) 


The release of solute into the liquid during 
solidification can be described by considering an 
average concentration in an arbitrary control volume 
which is undergoing phase change (Voller, Brent and 
Prakash 4 ). This control volume can be treated as 
partially solidified with an average concentration 
defined as: 


C = /,C, + / i C ( 


(17) 


Using (12), the energy equation (3) can be written: 

P[C‘(T)^+cV.(vt)j = XV 2 T. (13) 
at 9 


where f, = 1 - /, is the local solid volume fraction. 
Since diffusion in the solid is neglected, the 
concentration in the solid remains constant over time. 
Noting that C, = kC, we can thus write: 


To solve equation (13), an effective heat capacity 

coefficient 8 /, /dT has to be calculated. We define — =—-^-(1—11: )C, +( 1 — /, )— — J - (18) 

dt dt dt 


df t df t /dn( /, )J„ J /, ),T X H /, ),T, 

3T ST/dn T 2 T 2 +T 2 

where the subscripts n (denoting the normal direction), 
x and y denote differentiation. 

Since isothermal phase change is under consideration, 
the liquid fraction undergoes a step change when the 
interface crosses a grid line. This abrupt change in the 
liquid fraction, defined by the step function (10), can 
cause serious numerical instabilities. To overcome 
this problem, a control volume was defined around 
each grid point, in which the liquid fraction could be 
estimated. Phase change was considered to take place 
over one control volume, in which the step function 
(10) is replaced by a linear approximation: 

for TV > T m +AT /,=! 

T — T + AT 

f° r T m - AT < T/j <T m + AT /, = u ^ (15) 

for T,j <T m - AT /, = 0 


When (18) is used in the solute transport equation (4), 
we obtain the solute conservation equation in the 
form: 

^-+V.( VC i )=DV 2 C i +S c (19) 

dt 

in which 

S c ^(l-k)C l+ f,^- ( 20 ) 

The formulation for solute transport during phase 
change described by equations (19) and (20) allows 
for the solution for liquid concentration only and 
hence bypasses the concentration discontinuity at the 
interface. 

NUMERICAL METHOD 

An algorithm entitled SOLCON f 6 , which incorporates 
the closely coupled solution of the transport equations 


f SOLidification and CONvection 
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in the vorticity-stream function formulation, was used. 
To ensure stability of the computational process, all 
source terms and non-linear coefficients depending on 
liquid fraction are linearized based on the value of 
liquid fraction obtained from the previous iteration. 

The vorticity, stream function and energy equations 
were discretized using central differences and solved 
by a modified ADI scheme with internal iterations. 
Interface boundary conditions for vorticity and stream 
function were applied at those mesh points in the solid 
sub-region which are adjacent to the liquid. For the 
calculation of vorticity boundary conditions, the 
definition of vorticity was used: £ = VxV. The 

boundary condition y/ = 0 was used for the stream 
function. The concentration equation (19) was 
discretized and solved using a control volume 
approach. This ensures mass balance during phase 
change in the partially solidified control volume. A 
second order upwind scheme (SOU) was used for the 
convection fluxes with central differences for the 
diffusion terms. 

To account for the fact that the computed 
concentration is a cell average value, an exponential 
extrapolation procedure based on the liquid fraction 
has been introduced to find the values of the 
concentration at the solid/liquid interface. The liquid 
side interface solute concentration Q can then be used 
to determine the concentration in the solid as it forms. 

CODE VALIDATION 


To validate the code computations were performed for 
directional solidification of gallium-doped germanium 
crystal grown by the Bridgman-Stockbarger technique 
as described by Alexander et al. 2 . Comparisons were 
made for the case of sinusoidal acceleration with an 
amplitude of 10‘ 3 g and a frequency of 10' 1 Hz oriented 
parallel to the solid-liquid interface. 

In Alexander et al. 2 a pseudo-steady model was 
adopted with the assumption that the ampoule 
translation rate and the growth rate are equal. In their 
model the solid-liquid interface is located at a fixed 
distance from the top of computational domain, which 
is completely occupied by the melt. The aspect ratio 
of the computational domain was equal to 1. On the- 
other-hand the model used in the current study 
considers transient effects and hence the gradual 
decrease in the length of melting zone. Both solid and 
liquid phases are included in the computational 
domain (see Figure 1). 

The temperature profile is translated with a constant 
pulling velocity along the boundary causing the 
interface movement inside the domain. The boundary 
temperature profile, size of the computational domain 



Figure 1 . Model used in the present work for 
comparison with Alexander et al. 2 


and physical properties of the alloy were chosen to 
approximate the idealized model of the Bridgman- 
Stockbarger system in Alexander et al. 2 . T c was equal 
to 1 131 K and T h was 133 IK. The length of the 
computational domain was taken to be 21 mm, the 
height was 10 mm and the adiabatic zone was 5 mm. 
A uniform square 51 x 106 mesh was used. 

A steady solution for the temperature and flow was 
used as the initial condition for the transient growth. 
This steady solution was obtained by keeping the 
boundary temperature profile stationary with the solid- 
liquid interface located at 7 mm from the left wall of 
the ampoule. The initial solute concentration in the 
liquid was uniform at 1 at%. Solidification was first 
started with a constant gravity level of 1 jig and when 
3 mm of material was solidified, a sinusoidal 
acceleration with an amplitude of 10' 3 g and frequency 
of 10’ 1 Hz (oriented parallel to the solid-liquid 
interface) was imposed for another 2 mm of 
solidification. The translation velocity was 6.5 pm/s. 



Figure 2. Maximum vertical velocity for a 
sinusoidal acceleration with an amplitude 
of 10* 3 g and frequency of 10 1 Hz. 
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Figure 2 shows the maximum vertical velocity as a 
function of time computed using the present model. 
The velocity field is in phase with the residual 
acceleration. This result is in excellent agreement 
with Figure 7b from Alexander et al. . 


Table 1 . Comparison of the two models. 



U m „ 

(mm/s) 

v 

v max 

(mm/s) 

Yc n>«" Yo 
(%) 

Alexander etal. 2 

0.1140 

0.1325 

19.14 

Present work 

0.1175 

0.1338 

19.80 

% difference 

3.07 

0.98 

3.5 


Table 1 shows a quantitative comparison of the results 
obtained from the two models. U max and V max are the 
maximum horizontal and vertical components of 
velocity vector (during one period), and Yc is the 
segregation defined by 

Y< = (C M - CV )/ x 100% (21) 

where the three values of the concentration are taken 
at the interface and y 0 is segregation before the g-jitter 
starts. In the case of Alexander et al 2 y 0 was equal to 
zero, in the present calculations the initial segregation 
was 4.0%. 

The maximum change in segregation during g-jitter 
was equal to 19.14% (Alexander et al. 2 ) and 19.8% 
(our calculations). This was reached after 230 seconds 
of solidification. 

It is obvious that computed results are in very good 
agreement despite the difference in the physical and 
mathematical models. 


SOLIDIFICATION OF BI-SN ALLOY 

Simulations were performed for directional 
solidification of Bi-lat% Sn alloy in a Bridgman 
furnace. Property values are taken from Timchenko, 
Chen, de Vahl Davis and Leonardi 6 . 

The ampoule studied has a height of 6 mm and a 
length of 42 mm. The boundary temperature profile 
imposed on the liquid boundary consisted of a cold 
zone ( T c = 50 °C), linear temperature profile with a 
gradient 20 K/mm (for a length of 32.5 mm) and a hot 


zone ( T h = 700 °C). That is, conduction in the 
ampoule wall was not considered. The computational 
domain initially contains only liquid with a uniform 
solute concentration C 0 of 1 at% and uniform 
temperature of 700 °C. At the left boundary an initial 
temperature of 272 °C was imposed. The pulling 
velocity - the rate of translation of the boundary 
temperature distribution - was 3.34 |im/s, and 
solidification occurred from left to right as time 
progressed. 

Mesh validation 

To ensure the accuracy of the solution a mesh 
validation has been performed for 50 seconds of 
solidification. The amplitude of the gravity vector was 
taken to be 1 |ig, acting in a direction normal to the 
axis of the ampoule. 


Table 2. Mesh validation 



31 x 21 1 

31 x 421 

61 x 421 

mm/s 

2.286x10'' 

(0.7%) 

2.286x10"* 

(0.7%) 

2.302xl0 J 

(-) 

V 

mu 

mm/s 

1.395x10'' 

(0.35%) 

1.400x10"* 

(0.0%) 

1.400x1 O'* 
(-) 


Three different mesh sizes were used, with the number 
of mesh points equal to 31 x 211, 31 x 421 and 
61 x 421. The time step used with the 31 x 21 1 mesh 
was equal 0.1 s, and for the 31 x 421 and 61 x 421 
meshes it was chosen equal 0.01 s. The difference 
between those results is listed in Table 2. 

From Table 2, it can be concluded that a 31 x 211 
mesh can be used in the calculations. 

Solidification with g-iitter 

The effects of periodic disturbances on the 
compositional profile and segregation at the interface 
have been investigated. The amplitudes were varied 
from 10* 5 g to 10 2 g for a range of frequencies from 
10' 2 to 1 Hz. 

At first, computations were performed for 
solidification with a constant gravity level of 1 jig. 
After 1500 seconds, sinusoidal accelerations oriented 
parallel to the solid-liquid interface were imposed and 
computations with g-jitter were performed for a 
further 500 seconds. 


5 





Table 3. Summary of results 


Amplitude 

(m/s 2 ) 

Frequency 

(Hz) 

Maximum velocity (mm /s) 

Maximum segregation 

(OL\ 

Segregation at 500s 
(%) 


v m „ 

{/o) 

log 

0.01 

2.50 

0.89 

188.7 

90.0 

0.05 

1.26 

0.78 

192.4 

188.0 

0.1 

0.75 

0.32 

39.9 

39.9 

0.5 

0.19 

9.66x10 2 

20.2 

20.2 

1.0 

9.78x10° 

4.81x10° 

10.9 

10.9 

10°g 

0.01 

0.22 

0.11 

38.2 

33.0 

0.05 

0.12 

6.75x10° 

11.3 

10.3 

0.1 

7.49x10 2 

4.13x10° 

8.9 

8.9 

0.5 

1.90x10 2 

9.78x10° 

3.1 

3.1 

1.0 

9.95x10’ 

4.94x10° 

2.5 

2.5 

10°g 

0.01 

2.15x10° 

1.15x10° 

4.9 

4.3 

0.05 

1.27x10 2 

6.86x10° 

2.7 

2.7 

0.1 

7.68x10° 

4.23x10° 

2.3 

2.3 

0.5 

2.06x10° 

1.08x10° 

1.8 

1.8 

1.0 

1.14x10° 

6.02x10° 

1.8 

1.8 

10 g 

0.01 

2.41x10° 

1.27x10° 

2.0 

2.0 

0.05 

1.47x10° 

7.96x10° 

1.8 

1.8 

0.1 

9.66x10° 

5.29x10° 

1.8 

1.8 

0.5 

3. 85x10^ 

2.12x10^ 

1.8 

1.8 

1.0 

2.97x10° 

1.65x10° 

1.8 

1.8 


steady 

2.28x10° 

1.20x10^ 

1.8 

1.8 


For comparison purposes, computations with a steady 
gravity level of 1 jig were performed. The results of 
these computations are summarised in Table 3. It can 
be seen that for a given acceleration magnitude, larger 
segregation at the interface occurs for smaller 
frequencies or longer periods of disturbances. The 
same effect was reported in Alexander et al\ 
Disturbances with amplitudes of 10' 5 g produce very 
little effect on the segregation or the compositional 
profile. Disturbances with amplitudes of lO^g cause 
increases in the segregation between 2.7 to 4.9% 
(1.8% for steady lpg gravity) when the frequencies of 
the disturbances become lower than 0.1 Hz. Larger 
effects on segregation were observed for disturbances 
with an amplitude of 10 3 g. 

Figure 3 shows the vertical velocity at a reference 
point as a function of time, for frequencies of the 
disturbances of (a) 0.01 Hz, (b) 0.1 Hz and (c) 1 Hz 
with an amplitude of 10‘ 3 g. The reference point is 
moving with the interface and is always located 2 mm 


in front of the interface at the mid-height of the 
ampoule. The development of the velocity field occurs 
in phase with the gravitational acceleration. The 
maximum velocity exhibits a transient before reaching 
steady oscillations. This transient includes number of 
periods which increases with increasing frequency as 
can be seen clearly in Figure 3(c). As the frequency is 
reduced the maximum velocity increases. 

Figure 4 shows the velocity field for one cycle of the 
disturbance with an amplitude of 10' 3 g and frequency 
of 0.01 Hz, starting from 300 seconds. During the 
first half of the cycle the gravity vector points 
downwards, with the maximum amplitude occurring at 
325 seconds. Near the interface the predominant flow 
direction is also downwards. The velocity reaches a 
maximum at 330 seconds and then decreases as the 
magnitude of the gravitational acceleration decreases. 
At 350 seconds gravity reverses and as a result the 
velocity field also reverses, reaching the next 
maximum in about 380 seconds. 
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(c) 

Figure 3. Vertical velocity at the reference point 
for an amplitude of 10' 3 g and frequencies of 
(a) 10‘ 2 Hz, (b) 10 1 Hz and (c) 1 Hz. 


t - 325 tcc 



Figure 4. Velocity Field: one cycle of the 
disturbance with an amplitude of 10’ 3 g 
and frequency of 0.01 Hz. 


t^vrl I $ 3 7 9 II t$ 13 

I.Oux C: DOS 1.04 2 OS l.OJ 4.02 5 01 6 01 7.00 



L^vd I 3 5 7 9 H 13 15 

( - 100 tec C 003 I.IO 1 IS 3.20 4.25 S 30 6.35 7 40 



Lcvd 1 3 5 7 9 1 1 1 3 1 3 

I -300 toe Cl 0.03 I IB 2 12 3-45 4.31 5.72 « IS 7 99 



Level I 3 3 7 9 tl 13 13 

t - 300 wc C: 0.03 1.27 2.49 3.71 4.93 6.15 7.37 I 31 



Figure 5. Solute Field for the disturbance with an 
amplitude of 10 3 g and frequency of 0.1 Hz. 

Concentration contours at 0, 100, 300 and 500 seconds 
after the start of g-jitter, with an amplitude of 10 3 g 
and frequencies of 0.1 and 0.01 Hz are shown in 
Figures 5 and 6 respectively. For a frequency of 
0.1 Hz convection is weak and the concentration 
contours shown in Figure 5 have only slightly 
distorted. The segregation in this case increases from 
1.8% to 8.9% by the end of 500 seconds of 
solidification. For 0.01 Hz the maximum segregation 
reaches 38.2% after 188 seconds of solidification and 
then decreases to 33% at the end of 500 seconds. In 
this case the flow carries low concentration solute 
from the bulk of the liquid to the interface and then 
redistributes it such that high concentration solute is 
located in the upper and lower parts of the interface. 
The concentration contours become much more 
distorted compared to the higher frequency case. 
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L«*«l I 3 5 7 9 II IS 15 

C: 0.03 1.04 2.03 3.03 4.02 3.01 6.01 7.00 


L«v«l I 3 S 7 9 11 13 13 

I» 100 Me C: 0 03 1.16 2.27 3.39 4 50 5.61 6.73 7 44 



Figure 6. Solute field for the disturbance with an 
amplitude of 10 3 g and frequency of 0.01 Hz. 



Figure 7. Distribution of solute concentration 
at the mid-height of the ampoule after 
500 seconds of solidification. 


Figure 7 shows the solute concentration distribution at 
the mid-height of the ampoule at the start of g-jitter 
and after a further 500 seconds of solidification for the 
three frequencies of 0.01, 0.1 and 1 Hz, for 
disturbances with an amplitude of 10‘ 3 g and also for a 
steady gravitational acceleration. 

By the time g-jitter starts the first 5 mm of the sample 
had been solidified, creating a solute rich boundary 
layer in front of the interface. This decays 
exponentially to Cq. The peak value of concentration 
at the interface caused by solute rejection into the 
liquid reached almost 7 at%. 



Figure 8. Average concentration at the interface 
for the disturbances with an amplitude of 10 3 g 
and frequencies of 10 2 , 10' 1 and 1 Hz. 



Figure 9. Segregation at the interface for the 
disturbances with an amplitude of 10 3 g and 
frequencies of 10' 2 , 10 1 and 1 Hz. 


After a the further 500 seconds of solidification the 
maximum concentration at the mid-height of the 
ampoule reached 8.33 at% for the steady gravitational 
acceleration, 8.43 at% and 8.38 at% for frequencies of 
0.1 and 1Hz, and 7.26 at% for 0.01 Hz. The first three 
values only vary by 1.2% and correspond to diffusion- 
controlled plane front solidification. The largest 
difference in the maximum value (12.8%) occurs at 
the lowest frequency of 0.01 Hz. In this case 
convection has developed, resulting in the 
redistribution of solute in the melt and hence at the 
interface. The difference between the low 0.01 Hz 
frequency and the other cases can clearly be observed 
in Figure 8 which shows the history of the average 
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concentration at the interface during 500 seconds of 
solidification. 

The resulting segregation at the interface can be seen 
in Figure 9. Clearly the segregation at the interface 
fluctuates with a frequency corresponding to that of 
the disturbances. For a frequency of 0.1 Hz, the 
segregation exhibits periodic oscillations after 200 
seconds of g-jitter, whereas for a frequency of 0.01 Hz 
the periodic regime has not been reached by 500 
seconds. The different magnitude of the peaks in the 
case of the 0.1 Hz frequency and their periodic change 
are not real but rather an artefact of the fixed grid 
finite volume formulation which is used. In this 
formulation all the computed values are cell averaged 
values which are changing while the interface passes 
through one partially solidified cell. 

Finally, disturbances with an amplitude of 10 2 g are 
considered. This case exhibits not only quantitative 
change in the segregation for different frequencies but 
also different regimes of solidification. The time 
history of the average concentration at the interface 
during 500 seconds of solidification is shown in 
Figure 10. The response of the average at the 
interface concentration changes character with 
decreasing of frequency. At 1 Hz diffusion controlled 
growth is observed, with a gradual concentration 
increase due to release of solute in the liquid. 
Segregation is also increasing, reaching 10.9% after 
500 seconds of solidification. 

At 0.1 Hz the pattern is changing. The flow at first 
carries low concentration solute from the bulk of the 
liquid closer to the interface, some mixing occurs in 
the liquid and as a result the average concentration at 
the interface decreases while the interface is moving 



Figure 10. Average concentration at the interface 
for the disturbances with an amplitude of 10' 2 g 
and frequencies of 10 2 , 10 1 and 1 Hz. 


L«v*l 1 3 S 7 9 11 13 IS 

C: 0.05 1.27 2.49 3.71 4 93 6 15 7.37 6 59 



(•) 

L*v*l 1 3 5 7 9 11 13 15 

C 0.05 0.99 1.92 2 66 3.80 4.74 5.66 6.62 



(b) 

L*v*l 1 3 5 7 9 11 13 15 

C: 0.05 0.42 0.80 1.16 1.56 1 94 2 32 2 70 



(C) 


Figure 11. Solute concentration at 500 seconds 
for the disturbances with an amplitude of 10 2 g and 
frequency of (a) 1 Hz, (b) 0.1 Hz and (c) 0.01 Hz. 


through the liquid. However, once the flow starts to 
oscillate backwards and forward, the average 
concentration reaches a steady state. The segregation 
in this case reaches 39.9% by 500 seconds. At 0.01 
Hz the average concentration drops quickly due to 
mixing with solute from the bulk of the cavity. At 60 
seconds the flow reverses and brings back high 
concentration solute causing an increase in the average 
concentration at the interface. As the magnitude of 
the velocity increases in the opposite direction, further 
mixing occurs in the liquid and the average 
concentration drops to 1.5%. 

Figure 1 1 shows contours of concentration at 500 
seconds for the disturbances with an amplitude of 
10 2 g and frequency of (a) 1 Hz, (b) 0.1 Hz and (c) 
0.01 Hz. 

Figures 12 and 13 show the transformation of the flow 
structure from one configuration to another as a result 
of gravity reversal. These transitional stages are 
shown for 10" 3 g in Figure 12 and for 10‘ 2 g in Figure 
13; the frequency of the disturbance was 0.01 Hz. 
Times shown are the 52 nd to 56 th seconds of a 100 
second period. Gravitational vector reversal occurred 
at the 50 th second. 

For 10 3 g (Figure 12), the flow goes through a 
transition from a counter clockwise main cell to the 
development of two small secondary and one stronger 
main cells and finally a single clockwise cell. 

In the case of 10 2 g (Figure 13), distortion of the main 
flow cell has started after 30 seconds at which time the 
magnitude of the acceleration has started to decrease. 
The flow structure goes from a counter clockwise cell 
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structure to a flow with three secondary cells, then to a 
two cells structure and finally to the single cell flow 
structure at 57 th second. 


— 0.02 mm/s 



Figure 12. Flow reversal for the disturbance with an 
amplitude of 10' 3 g and frequency of 0.01 Hz. 
Times from top to bottom are the 52 nd to 56 th 
seconds of a 100 second period. 


0. 2 mm/s 
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Figure 13. Flow reversal for the disturbance with an 


amplitude of 10‘ 2 g and frequency of 0.01 Hz. 
Times from top to bottom are the 52 nd to 56 th 
seconds of a 100 second period. 


CONCLUSIONS 

We have investigated numerically effects of sinusoidal 
disturbances with amplitudes from 10' 5 g to 10' 2 g and 
frequencies from 10' 2 to 1 Hz on the solute 
redistribution and segregation at the interface during 
directional solidification of Bi-lat% Sn alloy. 

It was found that for large frequencies a higher 
amplitude of the gravitational acceleration is required 
to produce an effect on the segregation. For example, 
disturbances with frequencies from 0.5 to 1 Hz and 
amplitudes less then 10' 2 g produce very little effect on 
the segregation. For this frequency range an 
amplitude of 10‘ 2 g resulted in 20.2% and 10.9% 
segregation compared with a segregation of 1.8% for 
the steady 1 pg case. For frequencies from 0.05 Hz to 

0.1 Hz, an amplitude of 10' 3 g results in the segregation 
changing to 1 1.3% and 8.9% respectively. The largest 
effect on the segregation was produced by 
disturbances with a frequency of 0.01 Hz, where the 
maximum segregation was equal to 4.9% for an 
amplitude of 10\ 38.2% for an amplitude of 10 3 g 
and 188.7% for an amplitude of 10' 2 g. In the last case 
complete mixing of the solute in the cavity was 
observed. 
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